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We re-examine the well-known Wong formula for heavy-ion fusion cross sections. Although this 
celebrated formula yields almost exact results for single-channel calculations for relatively heavy 
systems such as ^®0-|-^‘^'^Sm, it tends to overestimate the cross section for light systems such as 
We generalise the formula to take account of the energy dependence of the barrier param¬ 
eters and show that the energy-dependent version gives results practically indistinguishable from 
a full quantal calculation. We then examine the deviations arising from the discrete nature of the 
intervening angular momenta, whose effect can lead to an oscillatory contribution to the excita¬ 
tion function. We recall some compact, analytic expressions for these oscillations, and highlight 
the important physical parameters that give rise to them. Oscillations in symmetric systems are 
discussed, as are systems where the target and projectile identities can be exchanged via a strong 
transfer channel. 
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I. INTRODUCTION 

It has been known for many years that fusion cross 
sections for several light heavy-ion ^sterns, such as 

i2c-pi2c (Ref. Q), (Ref. ( 3 ) and 

(Ref. 0), exhibit an oscillatory structure as a function of 
the incident energy. Two recent papers by Esbensen Q 
and Wong [sj interpret the oscillations as due to the ad¬ 
dition of successive individual partial waves as the energy 
increases. The effect is of course most important for iden¬ 
tical spin-zero nuclei, because in that case the odd partial 
waves are totally absent and the relevant energy spac¬ 
ing between successive contributing angular momenta is 
consequently much larger. Both of these authors discuss 
this phenomenon in the context of the well-known Wong 
fusion cross section @ derived from the Hill-Wheeler ex¬ 
pression 0 for the penetration of a parabolic potential 
barrier. See also Ref. @ for a recent publication in a 
similar context. 

In fact the above interpretation was first proposed 
some 30 years ago by Poffe, Rowley and Lindsay @ who, 
furthermore, gave a compact and accurate analytic ex¬ 
pression for the oscillations that displays succinctly the 
dependence on the relevant physical parameters of the 
system. Our purpose in this paper is to discuss the 
derivation of the Wong cross section - both its smooth 
and oscillatory terms - and to present some inadequa¬ 
cies of the standard fusion formula that treats the bar¬ 
rier height, position and curvature [R, i?B,/iw] as inde¬ 
pendent of the incident energy E and thus, implicitly, 
independent of the angular momentum 1. 

This re-analysis of the qualities and weaknesses of the 
Wong formula is important in order to be able to distin¬ 
guish between discrepancies arising from the approxima¬ 
tions used in its derivation and those arising from addi¬ 
tional physical effects such as entrance-channel couplings 
or limitations on compound-nucleus formation. 

The paper is organized as follows. In Sec. II, we derive 


the Wong cross section and propose an extended formula 
which takes into account the energy dependence of the 
barrier parameters. We show that the effect of the energy 
dependence is significant, particularly for light heavy-ion 
systems, and that the generalised formula reproduces well 
a full quantal calculation. In Sec. Ill, we discuss the os¬ 
cillatory part of fusion cross sections. We re-analyze the 
compact formula for these oscillation using the energy- 
dependent version of the Wong formula. We discuss the 
fusion of two identical nuclei as well as fusion between 
similar nuclei. Fusion oscillations in heavier systems are 
also discussed. In Sec. IV, we present our analyses of the 
experimental data for the and sys¬ 

tems. We summarize the paper in Sec. V. 


II. WONG FORMULA FOR FUSION CROSS 
SECTIONS AND ITS IMPROVEMENT 

In this section, we compare the results of the Wong 
expression for the fusion cross section that contains the 
three parameters [R, (see below) and the fusion 

cross section coming from a quantal calculation with a 
specified real potential and an absorption that is essen¬ 
tially black box. This can be achieved either with an 
ingoing-wave boundary condition or with an appropri¬ 
ately chosen imaginary potential. Here we choose the lat¬ 
ter method, checking carefully that our results are insen¬ 
sitive to changes in the imaginary potential used. With 
a potential model, the barrier position Rb, its height R 
and its curvature huj are all fixed by the potential param¬ 
eters. Indeed it is a very general result that the nuclear 
potential is essentially exponential in the tail, the region 
in which the Coulomb barrier occurs, at least for rela¬ 
tively light systems with low ZxZ^. In that case, all three 
of the above parameters are determined by two potential 
parameters, the depth Vq and the surface diffuseness a. 

We choose to fix these parameters in the following. 
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where Rb is the position of the barrier. Writing Vi<i(r) = 
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Furthermore, we then have a Coulomb plus nuclear bar¬ 
rier 


B = V{Rb) 
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and this quadratic equation for Rb gives the barrier po¬ 
sition as 

Rb = \rc + (5) 

where Rc is just ZxZ^^ jB [i^. We see that once the 
barrier height B is given, the position depends only on 
the surface diffuseness of the potential, and this also fixes 
the strength of the nuclear potential. This is the rea¬ 
son that it is generally sufficient to quote the parame¬ 
ters [i?, a] in our coupled-channels (CCFULL [ 13 )' calcula¬ 
tions. It is of course also sufficient for simpler optical- 
model/classical calculations of the fusion cross section. 
This procedure is useful because the parameter best de¬ 
termined by the fusion data is B, and furthermore there 
are various prescriptions for obtaining a good theoreti¬ 
cal value for this quantity; for example the Bass [HI and 
the Akyiiz-Winther [l^ potentials that give very similar 
B values over a wide range of heavy-ion systems, both 
yielding B = 6.1 MeV for the system mainly 

discussed in this paper. 

As noted above, this reduces the three independent 
parameters of the Wong expression to two, thereby pro¬ 
viding a much more rigorous constraint on the physics 
of the problem, as we shall see below. Notice that the 
exponential potential tends to be deep, and such deep 
potentials have been advocated in Refs. fl3l - [l5l| . 

If all of the flux crossing the Coulomb barrier fuses, 
then the fusion cross section is given by 


a = 


2mE 


1=0 


( 6 ) 


and using the Hill-Wheeler formula for the transmission 
Ti through a parabolic barrier, we have @ 
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where in the barrier region, the potential is taken as 

V{r) = B + - {r — RbY■ (The quantity huj 

is the quantum of energy corresponding to the inverted 
barrier and is generally referred to as the ‘barrier cur¬ 
vature’ [see Eq. (US].) Taking [B, as fixed, and 

replacing the summation by an integral, one obtains the 
very influential Wong formula for the fusion cross section 
for a single potential barrier 


hUJ n , , 

O' = ^^B In ( 1 + exp 




( 8 ) 


which yields, in the limit E — B ^ huj/2Tr, the classical 
result 


Ea = 7ri?| {E - B). 


(9) 


The first and second derivatives of this classical equation 
and its quantal version yield 
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where x = {271 /huj){B — E) and 9 is the Heaviside step 
function. These functions (especially the second) have 
been extremely important in the development of the no¬ 
tion of a fusion barrier distribution [^, where the 
fact that c?2 {Ea) / dE^ is strongly peaked near to the bar¬ 
rier is a crucial point. 

It is important to notice that higher above the bar¬ 
rier, one must question the approximations that lead to 
Eq. (I8|). Figure [T] shows that for higher angular momenta 
the barrier in the full potential (Coulomb -|- nuclear -|- 
centrifugal) occurs at a separation Re smaller than its 
I = 0 value, Rb- Now the grazing angular momentum is 
given by 


Igilg + 1 )^^ 
2mR^ 


= E-Ve, 


( 12 ) 


where Ve is the sum of the Coulomb and nuclear poten¬ 
tials at Re- Of course the curvature changes too, not 
only because the position of the barrier has shifted, but 
also because there is now a centrifugal contribution: 
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where the second derivatives are evaluated at r = Re- 
A better approximation for the fusion cross section is 
now just 
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FIG. 1: For an incident energy E above the Coulomb bar¬ 
rier B, the radius for the barrier with grazing angular mo¬ 
mentum Ig occurs at a separation Re < Rb- The sum of 
the Coulomb and nuclear potentials at this point is Ve < B. 
Furthermore, the curvature &j_b of the barrier in the total 
potential is larger than its value for I = 0. The curves are 
calculated for an exponential potential with a = 0.8 fm and a 
depth that yields B = 6.22 MeV for the system. For 

Z = 0 this potential has Rb ~ 7.44 fm and Juje = 2.52 MeV. 


FIG. 2: (Color online) (a) Fusion cross sections for the 
system obtained with the Wong formula with Z = 0 
values of [B, /ku] (dashed curve) and with their energy- 
dependent values (dot-dashed curve). The Bose symmetry 
of the identical spin-zero system is ignored, and both even 
and odd partial waves are summed up in the cross sections, 
(b) The energy-dependent Wong result agrees extremely well 
with a full quantal calculation shown by the solid line. How¬ 
ever, the latter shows very weak oscillations (« 1 mb) coming 
from the partial-wave sum even though all Z are summed. 


which is simpiy the standard Wong formuia but with 
an energy-dependence [B,Rb,uj] —>■ [Ve, Re,^e] derived 
for the grazing an gui ar momentum Ig at the energy in 
question. See Ref. for a simiiar extension to the cias- 
sicai formuia, Eq. (0. 

Of course we are stiii making an approximation here, 
because it is ciear from Fig. [T] that the barrier height, 
position and curvature depend expiicitiy on 1. Of course 
if one does not assume the same vaiues of the parameters 
for aii Z, then the integrai over I ieading to Eq. da can¬ 
not be performed anaiyticaiiy. But there is nothing to 
prevent one choosing a different parameter set (for aii Z) 
at each energy. Cieariy the set of vaiues for Ig is the 
best choice. One then stiii obtains the compact expres¬ 
sion (HI though we shouid stress that one now needs 
to derive numericaiiy the parameters [Ve, Re,>-^e] (see 
Appendix A). 

We show in Fig. [2] some caicuiations for the system 
that we shaii concentrate on in this paper. 
For the moment, they ignore the Bose symmetry of this 
identicai spin-zero system and sum over aii even and 
odd partiai waves. The dashed iine in Fig. [2] (a) shows 
the Wong cross section with [B,Rb,uj] fixed at aii en¬ 


ergies to their vaiues for I = 0. That is [B,Rb,uj] = 
[6.22 MeV, 7.42 fm, 2.52 MeV], generated by a potentiai 
with a = 0.8 fm that gives this barrier height. The resuits 
are seen to be very different from the dot-dashed iine that 
uses the energy-dependent vaiues of these parameters. In 
Fig- El (b) we compare the iatter resuits with a quantum 
mechanicai caicuiation of the cross section. The scaie of 
the verticai axis has been chosen to emphasise the high 
quaiity of the fit when the parameters are energy depen¬ 
dent. However, it aiso shows that even with aii partiai 
waves, there are smaii osciiiations in the fuii quantai cai¬ 
cuiations. We shaii derive an expression for these in the 
next section. 

The foiiowing comments are appropriate at this point: 

• The approximation with parameters fixed at their 
Z = 0 vaiues is poor, especiaiiy for iight systems. 
For these systems, the Couiomb interaction is reia- 
tiveiy smaii and aiso the reduced mass m is smaii, 
therefore the centrifugai potentiai may piay a more 
important roie than the Couiomb potentiai. For 
heavier systems, on the other hand, the Couiomb 
potentiai is strong and thus the barrier is rigid 
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against a variation of angular momentum. In that 
situation, the conventional Wong formula is reason¬ 
able. Figure [3] shows the barrier position Rb and 
the barrier curvature fiw as a function of angular 
momentum I for two different systems, 
and ^^O-l-^^'^Sm. One can see that the variation is 
marginal for the heavier system, ^^O-l-^^'^Sm, while 
both Rb and hui change considerably as a function 
of I for the lighter system, 

• However, note that experimental fusion cross sec¬ 

tions (especially if they are not terribly precise) can 
frequently be fitted even for light heavy-ion sys¬ 
tems with parameters that are simply chosen to 
do so, and do not necessarily bear much relation 
to the physical potential. For example, we shall 
see in Sec. III-D that the experimental 

fusion cross section can be reasonably well fitted 
with [B, RByfiuj] = [5.6,6.3,3.0] (see Fig. [5] below). 

• The Fl-dependent Wong formula works well because 
the parabolic approximation is good, since the en¬ 
ergy E coincides with the barrier height for I = Ig. 
Even though this approximation breaks down at 
energies below the barrier, only a small number of 
partial waves in the vicinity of Ig contribute. The 
cross section is relatively large above the barrier 
and the resulting discrepancies are small. 

• At sub-barrier energies, there is no Ig as such, and 

the best one can do is to take the parameters for 
1=0, though one is already below the barrier for 
this value of 1. For higher I values, the situa¬ 
tion is worse, but taking an /-dependent value of 
the parameters does not improve matters, as the 
parabolic approximation is in any case intrinsically 
poor below the Coulomb barrier Here the 

cross section is small, so the errors are relatively 
more important. 


III. FUSION OSCILLATIONS 
A. Fusion of two identical nnclei 

In most situations, the value of hui is irrelevant at en¬ 
ergies above the Coulomb barrier, simply determining 
the rate of fall-off for E < B. (Note that in that re¬ 
gion the parabolic approximation itself will become in¬ 
adequate at sufficiently low energies [S[2l|.) However, 
(as first pointed out by Poffe et al. [9|, and more re¬ 
cently in Refs. 0; 0) there is a rather unique situation 
where huj is important above the barrier, principally for 
light symmetric spin-zero systems (though see below for 
other examples). Here, only even values of the grazing 
angular momentum are allowed and their barriers may 
be sufficiently well separated in energy for their succes¬ 
sive addition to the cross section to give rise to fusion 
oscillations. 



I 

FIG. 3: (Color online) The barrier position, Rb, (upper 
panel) and the barrier curvature, huj, (lower panel) as a func¬ 
tion of angular momentum, /, for the ^^C-|-^^C system (solid 
line) and the ^®0+^^^Sm system (dashed line). 

Let us now derive the earlier expression of Poffe et al. 
0 for these fusion oscillations. This is obtained by using 
the exact Poisson summation formula 

poo 

/ /(A)exp( 27 rTOa)dA. 

l m -^0 

(15) 

When applied to Eq. 0, the m = 0 term gives the usual 
Wong result da. The terms m = ±1, however, give 
rise to an oscillatory contribution from the poles of the 
transmission function above and below the real-A axis 
respectively. In the energy region E — B ^ huj/2'K, the 
nearest poles simply give (see Appendix B) 

Actosc = 2TTR%hujE exp(-2^) sin(27r/g), (16) 

where Ig is the continuous, energy-dependent grazing an¬ 
gular momentum of Eq. m, and the quantity ^ is given 

by 

IEje irmRl; _ tt ftiUE , . 

^~2lg + l' ~2dVEldlg' ^ ^ 

The resultant oscillations are shown on top of the smooth 
part of the cross section in Fig. [5] (b). They are of the 
order of 1 mb, and it is unlikely that any reasonable ex¬ 
periment would be able to observe these. Terms from 
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E (MeV) 

FIG. 4: (Color online) The symmetrised cross section (even I 
only) for is compared with the cross section for all 1. 

One sees an order of magnitnde difference in the oscillatory 
terms (see text). It is unlikely that the oscillations shown in 
Fig. 2 (b) (« 1 mb) could be observed in a system without 
some asymmetry between odd and even I values. 

|m| > 1 and from more distant poles merely introduce 
higher multiples of the exponent in the above equation 
and are clearly completely negligible. 

However, for a system of two identical spin-zero bosons 
(for example, ^^C-b^^C) the Bose symmetry forbids odd 
angular momenta, and the fusion cross section is given 
by twice the sum over the even partial waves. This can 
be achieved by including the factor 

1 -b [—1]* = 1 -b exp^iril) (18) 

in Eq. (USD. The 1 gives the usual non-symmetrised result 
but the extra term now gives an oscillatory contribution 
(from m = 0 and m = —1), 

Eaosc = 27rR^Hu!E exp(-^) sm(7rlg), (19) 

which is significantly larger than before, because the 
(negative) exponent has been reduced by a factor 2. This 
is demonstrated in Fig. 01 where the dashed curve rep¬ 
resents the sum of the smooth, energy-dependent Wong 
cross section and the above expression for the oscillatory 
part. The solid curve is the full quantal calculation, and 
it can be seen that the above expressions give an excel¬ 
lent approximation to the exact result. The dot-dashed 
curve shows the smooth part of the Wong cross section. 
One sees an order-of-magnitude increase in the oscilla¬ 
tions arising from symmetrisation. 

The above result is easily generalised to two identical 
spin-1/2 fermions such as ^^C-b^^C. Because the oscil¬ 
latory cross section for odd partial waves is clearly just 
minus that for the even ones, the total CTosc is simply 
given by the difference of the statistical weights for S' = 0 
(which gives even I with weight 1/4) and for S = 1 (giving 
odd I with weight 3/4). Thus we see that the oscillations 
will be reduced by a factor 2, and will have the opposite 


phase from the symmetric system. Similarly, the oscilla¬ 
tions for -b (1 0 1) will be reduced by a factor 
of 3 0. 


B. Fusion between similar nuclei: the role of elastic 
transfer 

The other interesting problem of this type is that of a 
system such as ^^C-b^®0, which has also been discussed 
earlier by Kabir, Kermode and Rowley [2^. Here, one 
can again consider the oscillations from even I and odd 
I separately. Without any further effect, they will cancel 
out. However, the presence of the elastic a-transfer chan¬ 
nel will introduce a parity dependence into the problem. 
This is most easily seen by considering the total elastic 
scattering, where to the amplitude fei{6) for direct elas¬ 
tic scattering of the we must add an exchange term 
/trans(7r — 6>), the amplitude for elastic transfer at the an¬ 
gle and energy in question. This yields a total amplitude 

/total(^) /el(0) + /trans(7r — 0 ), (20) 

and using P/(cos(7r — 6)) = (—l)*P;(cos0) we obtain 

Sf = S-f' -b (-l)'^/*™", (21) 

that is, different effective S'-matrix elements for the odd 
and even partial waves. This effect can be simulated (as 
in Refs. [23l - l28j ) by introducing a parity-dependent op¬ 
tical potential. However, to obtain a simple orientation 
of the effect of transfer (and again an analytic expression 
for the oscillations) we note that the ^-matrix for direct 
reactions is generally peaked around the grazing angular 
momentum and can be approximated by adS^^/dl. For 
relatively small a (which is of course limited by unitar- 
ity), this gives an effective S'-matrix 

S^«{l) = S^\l+[-lfa). (22) 

That is, the transfer dynamics results in a small shift of 
the original elastic S-matrix in opposite directions for odd 
and even I values. Clearly for a non-zero a the poles giv¬ 
ing rise to Eq. o come into play and indeed for a=l/2, 
we obtain the same magnitude oscillations with a shift in 
phase. The general result is that the trigonometric func¬ 
tion in Eq. dUD is replaced as follows 

^ . sin(7r[/g-ba])-sin(7r[Zg-a]) 

sin(7rig) -)> - - - 

= cos(7rZg)sin(7rQ;). (23) 

Of course in this expression a may be energy dependent, 
and indeed if one uses a parity-dependent potential, this 
will naturally give rise to a shift that depends on the 
grazing I and thus on E. Ref. demonstrated that it 
is possible to obtain a reasonable fit to both the large- 
angle elastic scattering and the fusion oscillations with 
the same parity-dependent potential. 
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C. Heavier systems 


The presence or absence of measurable fusion oscilla¬ 
tions arising from the symmetrization for identical sys¬ 
tems depends mainly on the quantity ^ in Eq. (1171) be¬ 
cause it appears in the exponent in Eq. dni). Using an 
exponential nuclear potential, Eq. m yields exactly 


hujE = 


2E - Ve - 


Re'- 


amRs 


(24) 


with £ = 4{E—Ve) + Vce, where Vce is the Coulomb po¬ 
tential at Re- Close to the s-wave barrier B, the nuclear 
potential is relatively small, and there e ~ AE — 3Ve- 
However, at higher energies and higher grazing I values, 
the barrier is pushed to lower radii and the nuclear po¬ 
tential gives the dominant contribution. Here Vqe can 
be neglected so that e ~ A{E — Ve) yielding 


TT [^/2E-Ve 4a 
2 \E-Ve 


(25) 


or in terms of B 

( (l- — 

^ ^ 2 V 2a V Re 


1 


2a E 
Re E — B 



Asymptotically {E ^ B) this reduces to 


(27) 


and this gives a good idea of the mass dependence of 
the magnitude of the resulting oscillations. For heavy 
systems, i?_E is larger, and so is Thus the oscillations 
will become more difficult to observe experimentally in 
heavier systems. 

The top panel of Fig. [5] compares the value of ^ for 
the ^^C-b^^C system (dot-dashed line) with that for 
^®Si-|-^®Si (solid line). For the former reaction, we have 
used the same potential as in Fig. [21 while we have used 
an exponential potential with [H, a] = [28.8 MeV, 0.8 fm] 
for the latter system. Even though the ^ parameter for 
^®Si-|-^®Si is larger than that for the ^^C-I-^^C system by 
a factor of only about 2, its effect on the amplitude of 
the fusion oscillations is dramatic, as seen in the middle 
panel of Fig. |5l The oscillations shown here are obtained 
from Eq. (HH) but they can also be simply obtained from 
the quantum mechanical sum over h because the oscilla- 
tons arising from symmetrization are exponentially larger 
that those with no symmetry, they are essentially given 
by the difference between these two cross sections. This 
reduces to 

^osc = ^ ( E (2^ + m - E (2^ + 1)^; j , (28) 

\ I even I odd / 


where we note that this result is true only if the sum over 
all partial waves is sufficiently smooth. 



E-B (MeV) 

FIG. 5: (Color online) (a) The ^ parameter defined by 
Eq. (I17II in the oscillatory part of fusion cross sections. The 
dot-dashed and the solid lines show the ^ parameters for the 
system the ^®Si-b^®Si system, respectively, as a func¬ 
tion of the energy measured from the barrier height B of each 
system. Both lines are obtained with an exponential potential 
with the diffuseness parameter of a=0.8 fm. (b) The oscilla¬ 
tory contribution of fusion cross sections given by Eq. (HU. 
(c) The total fusion cross sections for the two systems given 
as a sum of the smooth part and the oscillatory part. 


Of course experimentally it is not possible to sepa¬ 
rate the oscillations out of the total cross section, and 
although they stand out in the ^^C-b^^C system (see bot¬ 
tom panel of Fig. ED, they are not apparent in the total 
cross section for ^®Si-b^®Si. For this reason the relevant 
experimental data are sometimes represented in the form 
of d{E(T)/dE [J. This is a useful representation since the 
derivative of the smooth part of Ea is essentially con¬ 
stant. We shall look at this again below. 

In the middle panel of Fig. |5l we have seen that the os¬ 
cillations for the ^®Si-b^®Si system are largely damped out 
at energies less than around 10 MeV above the barrier, 
while they start much earlier, and are much stronger, 
for the ^^C-b^^C system. In order to understand this. 
Fig. [5] compares the derivative of the penetrability for 
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FIG. 6: (Color online) The first energy derivative of the 
penetrability for each angular momentum I, dTi/dE, for the 
system (upper panel) and for the ^®Si+^®Si system 
(lower panel). The numbers above the peaks are the corre¬ 
sponding values of 1. 


each angular momentum with respect to the energy, that 
is, dTi/dE, for the system (upper panel) with 

that for the ^®Si-|-^®Si system (lower panel). For the lat¬ 
ter system, the individual peaks are narrower (due to 
smaller values of feu), but they are pushed much closer 
together, making them less well resolved from the peaks 
for adjacent values of L In order to obtain the same 
resolution that one has for / = 12 in the ^^C-|-^^C sys¬ 
tem (that is, the same degree of overlap between adja¬ 
cent peaks), one has to go to an I value of around 36 
or bigger in the ^®Si-|-^®Si system. This corresponds to 
a much larger energy compared with the former system, 
and indeed the limiting angular momentum for fusion for 
^®Si-|-^®Si appears to be fe^ax = 38, according to Vineyard 
et al. [ 2 ^. 

Several other data sets for the ^®Si-|-^®Si system ex¬ 
ist I 3 O-I 33 and indeed one of the aims of Aguilera et 
al. |32l| was to search for oscillations in this system. Un¬ 
fortunately they were unable to conclude the existence 
of structures within the 3% statistical error in their ex¬ 
periment. Their data are represented in the form Ea 
in the upper panel of Fig. [7] (open circles), along with 
data from Ref. [s^ (solid circles). It can be seen that 
the normalizations of the two data sets disagree some¬ 
what, presumably simply due to systematic differences 
in the experiments. Also shown in the figure are two 
uncoupled optical-model calculations with a=0.6 fm and 
1.2 fm (both with 5=28.7 MeV). The slope of the former 


fits that of the Aguilera data but the larger diffuseness is 
required to fit the slope of the data of Ref. [1^. In the 
lower panel of Fig. [7] the two data sets and calculations 
are again shown but in terms of d{Ea)/dE. Neither data 
set has the precision to prove the existence of oscillations, 
but it is interesting to note the significant difference in 
the magnitude of the oscillations produced by the differ¬ 
ent values of the potential diffuseness a (see Eq. (EZD). 
Better data are clearly required to resolve this question, 
and to establish the existence or absence of any oscilla¬ 
tions due to symmetrization. 



E (MeV) 


FIG. 7: (Golor online) Two different sets of measurements [^ . 

for fusion in the ^®Si-l-^®Si system are shown as Ea in the 
upper panel and as d{Ea)/dE in the lower one. The data 
are compared with uncoupled calculations with two different 
values of surface diffuseness in the nuclear potential. See text. 

Very recently Stefanini et al. [ 3 ^ have made more de¬ 
tailed measurements for this system with a very small en¬ 
ergy step (around 0.25 MeV in the center-of-mass system) 
and with better statistics (better than 1%). Their ex¬ 
perimental data, when displayed as d{Ea)/dE, indicate 
distinct structures at energies only around 5 MeV above 
the barrier. They are much more pronounced than those 
shown for the uncoupled calculation with a=1.2 fm in 
Fig.m This observation is incompatible with the mech¬ 
anism discussed here, and the oscillations in the new 
experimental data must have a different origin. Their 
presence appears to be related to strong couplings to the 
low-lying collective modes that exist in the ^®Si nucleus, 
but also seems to depend on the nature of the absorptive 
potential used in the calculations iH. 

In order to look for the effects of channel cou- 


























plings, the ‘experimental barrier distribution ’ El, E3 
D{E) = (P{Ea)/dE'^ is frequently used, and in the 
present context one should remember that the function 
d{Ea)/dE can also have structures close to the unper¬ 
turbed Coulomb barrier due to couplings. The width of 
a typical barrier distribution is proportional to Z\Z 2 and 
will lead to obvious structures in d(Ea)jdE for heavier 
systems with strong collective modes. However, for heav¬ 
ier systems, the ‘symmetrization’ oscillations will not be 
measurable. The system ^®Si-|-^®Si may be a special 
intermediate system where both types of structure are 
present simultaneously. 


IV. COMPARISON WITH EXPERIMENTAL 
DATA ON LIGHTER SYSTEMS 

A. fusion reaction 

Let us now analyze the actual experimental data for 
fusion of carbon isotopes and discuss the observed fusion 
oscillations. We first discuss the fusion of the ^^C-|-^^C 
system using a single-channel optical model. 

We have already shown in Fig. [3] the energy- 
dependence of the parameters Re and Huje entering both 
the smooth part of the Wong cross section and its oscil¬ 
latory term for the potential of Fig. [T] that gives a good 
fit to the experimental -I- fusion cross section, 
including its oscillations. (Some of these energy varia¬ 
tions have also been discussed in Ref. E^, where they 
are expressed in terms of universal functions.) The func¬ 
tion dTi/dE shown in the upper panel of Fig. [S] is for an 
optical-model calculation with the same potential and 
various I values. Its width is 0.56 Huje E^, and this com¬ 
pares well with the values oihuiE from Fig. [31 taken at the 
appropriate peak energies (where the I values in question 
are grazing). The variation of all of these Wong parame¬ 
ters is very signihcant over the energy range in question, 
and it is important to remember that their I = 0 val¬ 
ues do not even fit the average data when inserted into 
Eq. ([T31) . 

In Fig. [51 the data are reasonably well fitted dill (at 
least up to around 20 MeV) with an energy independent 
set of [B,R,fiuj] = [5.6,6.3,3.0] but these parameters 
do not correspond to those coming from the potential- 
model fit shown in Fig. [3] (a). (Though the value of the 
oscillatory term is similar to the potential-model value 
in the region where the oscillations are important, be¬ 
cause the expression (O was used for it.) For instance, 
Eq. (|3|) indicates that this parameter set corresponds to 
an exponential potential with the diffuseness parameter 
of a = 2.0 fm, that is unphysically large. The Wong cross 
section can, therefore, be considered only as providing a 
simple parametrisation of the fusion data, although as 
pointed out earlier this is extremely useful if one remains 
close to the barrier. 

This point is well demonstrated by the failure to fit 
the cross section at energies above E « 25 MeV, and 



FIG. 8: (Color online) Parametrised ^^C-|-^^C cross sec¬ 

tion with [B, Rb, fiijj] = [5.6, 6.3, 3.0]. The smooth-only 
cross section is shown by the dashed curve. The experimental 
data are taken from Ref. [3^ . 




FIG. 9: (Color online) (a) Potential-model fit to the same data 
as in Fig.[8]with [B, a] = [6.22, 0.8] (even I only), (b) Same as 
(a) but with the transmission of Z = 14 reduced by a factor 2 
and all higher waves absent. 

Fig. m (b) shows that a likely explanation of this is the 
failure of higher partial waves to fuse. This fit is obtained 
with the transmission T;=i 4 reduced by a factor of 2, 
and all higher partial waves completely removed from 
the cross section. 

Such an assumption is not unreasonable, because in 
this region of the compound nucleus spin I, the excitation 
energy of the ^^Mg formed in the ^^C-|-^^C reaction at the 
barrier height corresponding to an angular momentum 
Z = /, is not sufficient for s-wave particle emission. That 
is, 

Q + Ebarrier(0 = i^*barrier(^ = 0 < i^yrast(/) + S,, (29) 

where Q is the reaction Q value (13.93 MeV), Eyrast 
is the yrast energy of ^"^Mg [i^, and Sx is the ap- 
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E (MeV) E (MeV) 

FIG. 10: (Color online) Fits to the cross section 

with different potential parameters. The surface diffuseness 
a is constrained both by the magnitude of the smooth cross 
section and by the presence of oscillations. The left panel 
(a) shows results with a — 0.6 fm and B — 6.54 MeV, while 
the right panel (b) has a = 0.9 fm and B = 6.06 MeV. An 
intermediate value of a = 0.8 fm appears best (see Fig. [S]). 


propriate particle separation energy (where x € neu¬ 
tron, S'„=16.53 MeV; proton, S'p=11.69 MeV; or al¬ 
pha particle, 5o,=9.31 MeV). But, for Z > 12 we find 
■^barrier = 0 > -^yrast (d - 1) -f S'a, permitting the emis¬ 
sion of an L = 1 alpha particle. Similarly for particle 
angular momentum L = 2, proton emission also becomes 
possible. However, such emissions are inhibited by the 
penetration of the corresponding centrifugal barriers as 
well as the relevant Coulomb barriers, and thus competi¬ 
tion with fission will become important in the spin range 
above /=12. 

In Fig.[T0]we show potential model fits to the data with 
different values of the surface diffuseness a. While a good 
fit was obtained with a = 0.8 fm as shown in Fig. [5] (b), 
Figs. [TO] (a) and (b) show the best hts with slightly 
smaller and larger values a = 0.6 fm and a = 0.9 fm, re¬ 
spectively. In each case, the value of B is adjusted slightly 
to obtain the best fit and a clear pattern emerges. If one 
decreases a this increases Rb (Eq. (O) which in turn 
increases the higher-E cross section. This can be com¬ 
pensated by a slight increase in the barrier height which 
correspondingly decreases the radius again. Thus there 
is a little ‘play’ in the value of the potential parameters 
but one cannot stray too far from the ‘best’ values with¬ 
out destroying the ht either at higher or lower energies. 
Furthermore, the presence of the oscillations provides an 
additional strong constraint. For the smaller a value of 
0.6 fm, the average cross section is htted very well, but 
the magnitude of the oscillations is strongly damped. For 
the larger a value of 0.9 fm, the magnitude of the oscil¬ 
lations is increased (possibly improved) but the smooth 
part of the cross section starts to deviate. Thus we see 
that the data provide strong physical constraints on the 
physical properties of the nucleus-nucleus potential that 



_I_I_ n _I_I_L_ 

0 10 20 30 

E (MeV) 


FIG. 11: (Golor online) Effect of channel coupling on the fu¬ 
sion excitation function for the ^'^G-I-^'^G system (top panel). 
The dashed line is the same as the solid line in Fig. [9] (b) ob¬ 
tained with the potential model. The solid line, on the other 
hand, is obtained by including the rotational coupling to the 
first 2"*’ state both in the projectile and the target nuclei. The 
depth parameter of nuclear potential is slightly adjusted in or¬ 
der to remove the trivial barrier renormalization due to the 
couplings. The penetrability Ti and its first derivative dTi/dE 
are shown in the middle and the bottom panels, respectively, 
for Z = 8,10 and 12. 


are lost when the parameters [B, Rb, huj] are regarded as 
independent variables. Thus while the elegance and sim¬ 
plicity of the Wong cross section should be recognised, so 
should its limitations. 

So far we have discussed the fusion oscillations for the 
system based on the single-channel potential- 
model calculations. In order to discuss the role of chan¬ 
nel couplings, the top panel of Fig. [TT] shows results of a 
coupled-channels calculation for the system. To 

this end, we include the rotational coupling to the first 
2 + state both in the projectile and the target nuclei with 
the deformation parameter of /3=-0.40 (with a radius 
parameter ro=1.06 fm). In order to remove the poten¬ 
tial renormalization due to the coupling [^ , we 
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slightly adjust the depth parameter of the nuclear poten¬ 
tial. We use a modihed version of CCFULL to solve the 
coupled-channels equations. The solid and the dashed 
lines in Fig. [TT] denote the results of the coupled-channels 
and the single-channel calculations, respectively. One 
can see that the main feature of the fusion oscillations, 
including the peak energies and the phase of the oscil¬ 
lations, are not affected much by the channel-coupling 
effects, and thus our discussions based on the simple po¬ 
tential model calculations remain valid. 

This conclusion is due to the fact that the barrier dis¬ 
tributions for non-zero I d, [H, [13 still show almost a 
single-peaked structure, as shown in Figs.fTTlfb') and (c). 
This originates for the following two reasons. Firstly 
the excitation energy of the 2^ state in is relatively 
large (i?2=4.44 MeV) and thus the adiabatic approxima¬ 
tion is good, l eading to the adiabatic-barrier renormal¬ 
ization do, [so, [^ ■ Secondly, with rotational coupling to 
an oblate nucleus, the lowest barrier, which is relevant to 
the adiabatic-barrier renormalization, carries most of the 
weight, as in the vibrational case. If had a prolate 
deformation, the lowest barrier would have the smallest 
weight in the barrier distribution, and the fusion oscil¬ 
lations would be much more affected by the couplings. 
In that case, the oscillations are signiHcantly damped, 
and thus the data essentially determine the sign of the 
deformation of 

B. fusion reaction 

Let us next discuss the ^^C-|-^^C systems, for which 
the effect of elastic neutron transfer is expected to play 
an important role. Based on the discussion in the previ¬ 
ous subsection, we shall use the potential model for our 
discussions. Fig.[T2](a) shows the fusion cross section for 
this system and its ht with the Wong formula (dashed 
line). Although the oscillation is less prominent than 
for the fusion cross section still shows signih- 

cant oscillations and the Wong formula does not account 
well for them. The solid line in Fig. [12] (a) shows a ht 
with Eq. (IT^ . Although the magnitude of the oscilla¬ 
tory structure is well reproduced, the oscillation is out of 
phase with the experimental data. Fig. |T2] (b) shows a 
ht with Eq. (l23t with 0 ;=—0.5. It is apparent that this 
ht is better than the other two, and thus the experimen¬ 
tal data appear to determine the relative phase of the 
transfer and elastic ^'-matrix elements. 

Naturally, the sign of a is related to the sign of a parity- 
dependent part of optical potential. In fact, if one intro¬ 
duces a parity-dependent barrier height, Ve ± (—1)*AE, 
it is easy to show that the grazing angular momentum Ig 
is given to hrst order in AV by 




E (MeV) 


E (MeV) 


FIG. 12: (Color online) (a) Comparison between the exper¬ 
imental fusion excitation function and theoretical curves for 
the ^^C-|-^®C system. The dashed line is obtained with the 
Wong formula, while the solid line shows a ht to the data us¬ 
ing Eq. dm). The experimental data are taken from Ref. [s^ - 
(b) Same as (a) but with Eq. (I23II and a=—0.5. 



FIG. 13: (Color online) Fusion excitation function for 
obtained with a parity-dependent potential. The 
parity dependence is introduced by replacing Vb with Vb(l -b 
( — l)*e) in the exponential potential ’Fjv(r) = where 

Vo takes a negative value. The solid and the dashed lines are 
obtained with e = —0.15 and 0, respectively. The angular 
momentum sum is truncated at 1 = 14 as in Fig. [9] (b). 


This implies that a is given by 

a = T(-l) 


/ 2mi?| 


(A - VeW ’ 


(31) 


Ig = =F (-1) 


AVe 


{E-VE)h^’ 


(30) 


where 


is the grazing angular momentum for AE=0. 


indicating that a is positive (negative) if the barrier 
height for even-partial waves is lower (higher) than that 
for odd-partial waves. The negative sign of a for the 
^^C-b^^C system suggests that the barrier for this sys¬ 
tem is lower for odd-partial waves. This is consistent 
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with earlier findings in Refs. In Fig.[T3 we show 

fusion cross sections obtained with a parity-dependent 
potential. The parity dependence is introduced by using 
an energy-dependent depth parameter, Vb(l + (“l)*e) in 
the exponential potential (here, Vq is defined as taking 
a negative value). The solid line in the figure indicates 
that the experimental data are well accounted for with 
e=—0.5, which indeed implies higher barriers for even 
partial waves. 

In contrast to ^^C-|-^^C, the quantity a has been found 
to be positive for the ^^C-|-^®0 system [2^, and we note 
that Baye [i^ has proposed a simple rule for the sign of 
the parity-dependent potential in terms of the mass num¬ 
ber ^core of the identical ‘cores’ of the colliding nuclei and 
the parities of the valence orbitals of the transferred 
nucleons. His rule is based upon a microscopic resonat¬ 
ing group method (RGM) within a two-center harmonic- 
oscillator shell model. Essentially the barrier for even 
partial waves is deemed to be higher (or lower) than that 
for odd partial waves if the quantity 

- (-1)"^“- TT, (32) 

i:valence 

is positive (or negative). In both of onr cases Hcore is 
even. For ^^C-|-^®0 the number of valence particles is also 
even and this expression is negative. For the ^^C-|-^^C 
reaction we have a single, odd-parity valence neutron 
(Pi/ 2 ), and the expression is positive. The sign of a de¬ 
termined from the fusion oscillations is consistent with 
these results for both of our systems. 

Note that the sign of the parity-dependent term, ob¬ 
tained from a fit to the angular distribution of elastic 
scattering, shows some ambiguity for several systems [dlj. 
Fusion oscillations may offer a direct and perhaps better 
way to determine the sign. We note, however, that even 
though the negative value of a fits well the fusion oscil¬ 
lations and is consistent with the Baye’s rule, a positive 
value appears more consistent with the data for energies 
around 13 MeV (see Fig. IT^ . It would be interesting 
to re-measure fusion cross sections in this energy region 
with higher precision in order to confirm whether there 
exists a shift in phase of the oscillations as a function of 
energy. 

V. SUMMARY 

The Wong formula has been widely used to estimate 
fusion cross sections for a given single-channel potential 
as well as to discuss the parameters which govern fu¬ 
sion. The formula is useful in discussing for example, the 
subbarrier enhancement of cross sections, providing ref¬ 
erence cross sections in the absence of channel conplings. 
For relatively heavy systems, such as ^^O-l-^^'^Sm, the 
formula indeed reproduces well the exact result except 
for the deep subbarrier region, where the parabolic ap¬ 
proximation itself breaks down. On the other hand, for 
light systems, such as the Wong formula tends 


to overestimate the cross section. In this paper, we have 
extended the Wong formula by including the energy de¬ 
pendence of the parameters entering the formula, that is, 
the barrier height, barrier position, and the barrier cur¬ 
vature. Evaluating these parameters at the grazing angu¬ 
lar momentum for each energy, rather than at Z = 0, we 
have shown that the energy-dependent version of Wong’s 
formula reproduces the exact result well, even for light 
systems. 

The symmetrisation of the system leads to a Wong 
cross section possessing an oscillatory contribution, for 
which a compact formula can be derived based on the 
parabolic approximation. We have shown that the for¬ 
mula for the oscillatory contribution can also be extended 
to the energy-dependent version. 

Fusion oscillations are most significant in light sym¬ 
metric systems with spin-zero nuclei. We have analyzed 
the experimental data for the ^^C-|-^^C system and have 
argued that the fusion oscillations provide a strong con¬ 
straint on the nuclear potential employed in a calcula¬ 
tion. We have also analysed the system, in 

which elastic neutron transfer again gives rise to oscilla¬ 
tions. We have shown that these oscillations are useful 
in determining the sign of the effective parity-dependent 
potential arising from elastic transfer. 

We have argued that fusion oscillations provide an im¬ 
portant tool for studying properties of the nuclear poten¬ 
tial, strongly-coupled channels at high excitation energy, 
and fission. This is especially true for the sys¬ 

tem, which plays an important role in several astrophys- 
ical phenomena and thus has been recognized as one of 
the key reactions in that domain. Understanding 

both the smooth and oscillatory parts of the cross sec¬ 
tion above the barrier will almost certainly be necessary 
in understanding it in the important astrophysical region 
well below the barrier. 
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Appendix A: Approximate formula for the 
energy-dependent barrier position 

For the exponential potential of Eq. the grazing 
angular momentum Ig is given by 

^ ^ a^i^2e^ ZxZ^e^ lg{lg + l)h^ 

Re 2mR% 

(Al) 

with 5R = Rb — Re, where Re is the position of the bar¬ 
rier for I = Ig (see Eq. (ini)). Because the first derivative 
of the total potential is zero at the barrier position, we 
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also have 

_ JR/a _ Z\Zie^ ^ lg{lg + l)h^ 

R\ rt-R% 

(A2) 


By combining these two equations, it can be shown that 
SR satisfies 


^SR/a ^ A - B + \Vcb[1 + SR/Rb - (1 - 5R/Rb)-^] 
B-IVcb{1 + SR/Rb) 

(A3) 

where Vqb is again the Coulomb potential at Rb, the 
position of the I = 0 barrier. 

Inserting SR = 0 on the right-hand side of this equa¬ 
tion, one obtains 


^SR/a ^ ^ 


E-B 
B - ^VcB 


(A4) 


and from this first-order approximation to SR we obtain 
(c.f. Ref. [13) 


Re = Rb - aln . (A5) 

Eq. (IA3I) can be easily iterated and rapidly converges to 
the exact result for Re- 

The top panel of Fig. [TT] shows a comparison between 
the exact value of Re (solid line) with those obtained 
in first order (Eq. (IASI) ) and second order (dotted and 
dashed lines respectively) for the system. The 

middle and the bottom panels show Ve and JEje used 
in the energy-dependent Wong formula (ITT)) , evaluated 
at the corresponding exact and approximate Re values. 
One can see that the first-order formula given by Eq. (IA5I) 
works well only at energies in the vicinity of the s-wave 
Coulomb barrier, where SR is small. The second-order 
result already leads to excellent results over the entire 
energy range shown, and higher orders coverge rapidly 
to the exact result. 


Appendix B: Derivation of oscillatory term 


In this appendix, we give a detailed derivation of 
Eqs. (IT^ and (1191) . First we note that the Hill-Wheeler 
formula for the transmission. 


T{E,X) = 


l + exp[|^(H-E;+S)]’ 


(Bl) 


has poles at the complex angular momenta A which sat¬ 
isfy 


exp 


27T 

hjj 


(^B-E + 


2mR2 ) 


= ^ (B2) 


with n = 0, ±1, ±2, • • ■. To find the poles nearest to the 
real axis with n = 0 and —1, we put Apoie = Sg + iXj, 



E (MeV) 

FIG. 14: (Color online) (a) exact and approximate barrier 
positions Re for the system as a function of the in¬ 

cident energy E. An exponential nuclear potential with the 
diffuseness parameter a=0.8 fm has been used. The solid line 
shows the exact value obtained with Eq. (IA3I) . while the dot¬ 
ted and the dashed lines are obtained from Eq. (lASl) in first- 
and second-order respectively, (b) as (a) but for the value of 
the s-wave potential at Re- (c) same for the curvature of the 
barrier for the grazing angular momentum Ig- 


and neglecting the second-order term in A/, we find from 
Eq. (IB2I1 that Xg and A/ satisfy 


and 


B-E + 


2mi?2 


= 0 , 


(B3) 


Xi 


hw mR? 
2 Xg h? 


(B4) 


In the vicinity of the poles, the denominator of the right 
hand side of Eq. m is given by 


1 -b exp 


27r 

tUjJ 


B - E + 

2mR2 ) 


27r Apole^ 
hu mR"^ 


(A Apoie); 


(B5) 


where we have used Eq. (IB2I) to derive this equation. 

We are now ready to evaluate the m = ±1 terms in the 
Poisson summation formula (1151) by contour integration. 
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From the contour enclosing the upper, right quadrant of 
the complex plane, we may write the to = 1 term as 


27r 

F 

27r 


ai = -— / Ar(F;, A)e^"^dA, 


{+) hu mR^ 277 iA'+,> 

^pole'^ 


(B6) 


0 


- / {iX)T{E,i\)e-^^^{idX) 


(B7) 


where Ap^j^ takes the positive sign in Eq. (lB4l) . Here, 
we neglect the second term of this equation; as T is real 
and < 1 on the imaginary axis, the integral is real and 
< (27r)“^ and merely gives a very small correction to the 
smooth part of the cross section. So one finds 




, 27rfia; niR^ 


g 2 iriAg 


— TT 

e 





(B8) 


Similarly for the to = — 1 term we have 


a-I 


. 2Trhuj mR? 


g- 27 riAg 



(B9) 


Combining Eqs. (1581) and (1591) and using sin(27rAg) = 
sin(27rlg + tt) = — sin(27rlg), one finally obtains Eq. (TTSl) . 

With identical particles the fusion cross sections must 
be symmetrized and are given by 


^ + (BIO) 

c\ m—oo „oo 

= p I] (-I)"/ AdAT(E;,A)e2™*^(lT*e”^), 

(Bll) 


where the positive sign in Eq. (IB10|) relates to the spa¬ 
tially symmetric case and the negative sign to the spa¬ 
tially anti-symmetric case. The symmetrization is seen 
to lead to two terms, to = 0 and to = — 1 , where the 
exponent is reduced by a factor 2. These terms now 
clearly dominate, and the fusion cross sections are ap¬ 
proximately given by 


27r r°° 

a = —j XdXT{E,\)ilTiE"^) 


± 


27r 


XdXT{E,X)e 


—inX 


(B12) 


One can evaluate these integrals as above, and using 
cos(7rAg) = — sin(7rlg), we obtain 


fJ — CTWong fJosc, 


(B13) 


where the considerably enhanced oscillatory term is now 
given by Eq. m- 

In the above proof, we have used energy-independent 
values for B, R and huj for simplicity of notation. How¬ 
ever, the results are easily generalised to the energy- 
dependent case, when the above total cross section is 
given by Eqs. da and da using Ve, Re and Huje- 
(More generally these quantities depend on the angu¬ 
lar momentum rather than the energy, see discussion in 
Sec. HI). 


[1] P. Sperr, T.H. Braid, Y. Eisen, D.G. Kovar, F.W. Prosser 
Jr., J.P. Schiffer, S.L. Tabor and S.E. Vigdor, Phys. Rev. 
Lett. 37, 321 (1976). 

[2] P. Sperr, S. Vigdor, Y. Eisen, W. Henning, D.G. Kovar, 
T.R. Ophel and B. Zeidman, Phys. Rev. Lett. 36, 405 
(1976); J.J. Kolata, R.M. Freeman, F. Haas, B. Heusch 
and A. Gallmann, Phys. Lett. 65B, 333 (1976). 

[3] B. Fernandez, C. Gaarde, J.S. Larsen, S. Pontoppidan 
and F. Videbaek, Nucl. Phys. A306, 259 (1978); J.J. 
Kolata, R.M. Freeman, F. Haas, B. Heusch and A. Gall¬ 
mann, Phys. Rev. G19, 2237 (1979). 

[4] H. Esbensen, Phys. Rev. G85, 064611 (2012). 

[5] C.Y. Wong Phys. Rev. G86, 064603 (2012). 

[6] C.Y. Wong, Phys. Rev. Lett. 31, 766 (1973). 

[7] D.L. Hill and J.A. Wheeler, Phys. Rev. 89, 1102 (1953). 

[8] C. Simenel, R. Keser, A.S. Umar and V.E. Oberacker, 
Phys. Rev. C88, 024617 (2013). 

[9] N. Poffe, N. Rowley and R. Lindsay, Nucl. Phys. A410, 
498 (1983). 

[10] K. Hagino, N. Rowley and A.T. Kruppa, Comp. Phys. 


Comm. 123, 143 (1999). 

[11] R. Bass, Nuclear Reactions with Heavy Ions (Springer- 
Verlag, New York, 1980). 

[12] O. Akyiiz and A. Winther, in Nuclear Structure and 
Heavy-Ion Collisions, Proceedings of the International 
School of Physics “Enrico Fermi, ’’Course LXXVH, 
Varenna, 1979, edited by R.A. Broglia et al. (North- Hol¬ 
land, Oxford, 1981). 

[13] N. Rowley, H. Doubre and C. Marty, Phys. Lett. 69B, 
147 (1977). 

[14] Y. Kondo, B.A. Robson and R. Smith, Phys. Lett. B227, 
310 (1989). 

[15] Y. Kondo, F. Michel and G. Reidemeister, Phys. Lett. 
B242, 340 (1990). 

[16] N. Rowley, G.R. Satchler and P.H. Stelson, Phys. Lett. 
B254, 25 (1991). 

[17] M. Dasgupta, D.J. Hinde, N. Rowley and A.M. Stefanini, 
Annu. Rev. Nucl. Part. Sci. 48, 401 (1998). 

[18] N. Rowley, A. Kabir and R. Lindsay, J. Phys. G 15, L269 
(1989). 










14 


[19] K. Hagino, N. Rowley and M. Dasgupta, Phys. Rev. C67, 
054603 (2003). 

[20] K. Hagino and N. Takigawa, Prog. Theo. Phys. 128, 1061 

( 2012 ). 

]211 N. Rowley, N. Grar and M. Trotta, Phys. Rev. C76, 
044612 (2007). 

[22] D.M. Brink, Semi-Classical Methods for Nucleus-Nueleus 
Scattering, (University Press, Cambridge, 1985). 

[23] A. Kabir, M.W. Kermode and N. Rowley, Nucl. 
Phys.A481, 94 (1988). 

[24] W. von Oertzen, Nucl. Phys. A148, 529 (1970). 

[25] W. von Oertzen and H.G. Bohlen, Phys. Rep. 19, 1 
(1975). 

[26] M. Lozano and A. Vitturi, Phys. Rev. C35, 367 (1987). 

[27] A. Vitturi and C.H. Dasso, Nucl. Phys. A458, 157 
(1986). 

[28] J.A. Christley, M.A. Nagarajan and A. Vitturi, Nucl. 
Phys. A591, 341 (1995). 

[29] M. F. Vineyard et al., Phys. Rev. C 41, 1005 (1990). 

[30] S. Gary and C. Volant, Phys. Rev. C 25, 1877 (1982). 

[31] Y. Nagashima et al., Phys. Rev. C 33, 176 (1986). 

[32] E.F. Aguilera et al., Phys. Rev. C 33, 1961 (1986). 

[33] A.M. Stefanini, G. Montagnoli, L. Corradi, S. Courtin, 
E. Fioretto, J. Grebosz, F. Haas, H.M. Jia, M. Mazzocco, 
C. Michelagnoli, T. Mijatovic, D. Montanari, C. Paras- 
candolo, F. Scarlassara, E. Strano, S. Szilner, D. Torres! 
and C.A. Ur, EPJ Web of Conf. 66, 03082 (2014). 


[34] I. Wiedenhover et al, Phys. Rev. Lett. 87, 14250 (2001). 

35 S. M. Harris, Phys. Rev. 138, B509 (1965); Phys. Rev. 
Lett. 13, 663 (1964). 

[36] N. Takigawa, K. Hagino, M. Abe and A.B. Balantekin, 
Phys. Rev. C49, 2630 (1994). 

[37] K. Hagino, N. Takigawa, M. Dasgupta, D.J. Hinde, and 
J.R. Leigh, Phys. Rev. Lett. 79, 2014 (1997). 

[38] P. Sperr et al., Phys. Rev. Lett. 37, 321 (1976). D.G. 
Kovar et al., Phys. Rev. C20, 1305 (1979). 

[39] M. Yasue et al., Nucl. Phys. A394, 29 (1983). 

[40] D. Baye, Nucl. Phys. A460, 581 (1986). 

[41] Y. Kondo, B.A. Robson and R. Smith, Nucl. Phys. A410, 
289 (1983). 

[42] T. Spillane et al., Phys. Rev. Lett. 98, 122501 (2007). 

[43] M. Notani et al., Phys. Rev. C85, 014607 (2012). 

[44] C.L. Jiang et al., Phys. Rev. Lett. 110, 072701 (2013). 

[45] For a <C Rc this gives the very simple relation Rb ~ 
Rc — a, though we shall not use this approximation in 
this paper. 

[46] The ground-state rotational band of ^“^Mg is well fit¬ 
ted up to the highest identified spin of Ref. M 
(/’^ = 10^ at 19.2 MeV) by the Harris parametrization ^ 
with 77o=2.22 ft^.MeV”^ and 77i=0.115 /i'*.MeV“®. This 
gives the energies E 12 +=‘25.80 MeV, E 14 +=33.05 MeV, 
i? 4 g+=40.85 MeV..., and so forth. 



